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Diamagnetism in wire medium metamaterials: theory and experiment 

I. Yagupov/ D. Filonov/ S. Kosulnikov/ M. Hasan/ I.V. lorsh/ and P.A. Belov^ 
^ITMO University, St. Petersburg 197101, Russia 

Strong diamagnetic response of wire medium with finite wire radius is reported. Contrary to the 
previous works where it was assumed that the wire medium exhibits only the electric response, we 
show that the non-zero magnetic susceptibility has to be taken into account for proper effective 
medium description of the wire medium. Analytical and numerical results are supported by the 
experimental measurements. 


I. INTRODUCTION 

Wire metamaterials are composed of arrays of opti¬ 
cally thin metallic rods embedded in a dielectric matrijsi. 
These structures have been extensively studied both 
theoretically and experimentally, due to their possi¬ 
ble applications in subwavelength image transferi^r.^ and 
control over spontaneous emission lifetime of quantum 
emitters^i^. While there are different geometries of the 
wire media, in our work we consider a special class, 
formed by a two-dimensional array of parallel wires of 
finite radius and length, shown schematically in Fig. [T] 
For the practical applications and specifically in order to 
model electromagnetic properties of wire metamaterial 
blocks of complex geometry, effective material parame¬ 
ters such as permittivity and permeability are required 
for proper homogenization. The effective dielectric per¬ 
mittivity of the array of parallel perfectly conducting 
wires has been obtained in a number of papersiSr— . Per¬ 
mittivity and permeability tensors are given byi^: 



FIG. 1. Geometry of the considered structure: array of per¬ 
fectly conducting wires of finite radius ro and length D is 
placed inside the dielectric matrix with permittivity Em- 


count for the diamagnetic response of these structures. 


II. THEORETICAL MODELLING 
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where I is the unity matrix, EmiRm are the permittivity 
and permeability of the host material, and e^z reads 
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First, let us consider the quasistatic approximation for 
the single unit cell. We assume that the uniform static 
magnetic field is aligned along x direction. Then, the 
surface currents originating at the surface of the wire 
are aligned along the wires and can be expressed as 
= —2fjLrnHsuY{(j)), where / is the azimuthal angle. 
Knowing the currents we can calculate the magnetic mo¬ 
ment of the unit cell M as J (i^r[r x j]. Dividing it by the 
unit cell volume will give us the magnetization density 
and the effective permeability. The integral can be taken 
analytically yielding: 


where ro is the wire radius, and d is the period of the 
structure. We note that the dielectric permittivity ten¬ 
sor is nonlocal, i.e. it depends on the wave vector kz- 
However, in the case of an infinite wire structure along 
the lateral direction, at normal incidence, we can consider 
only the tangential components of the dielectric permit¬ 
tivity which coincide with the dielectric permittivity of 
the matrix. Thus, at normal incidence, the light effec¬ 
tively interacts with uniform dielectric with permittivity 
Em and magnetic permeability equal to unity. However, 
it is clear that the surface currents at the wires interfaces 
should lead to nonzero magnetic susceptibility. In this 
work, we present analytical, numerical, and experimen¬ 
tal results to conclude that for the proper description of 
the wire media effective parameters it is required to ac¬ 
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To obtain the effective permittivity we use the results of 
the worbi^, where it is shown that wire medium belongs 
to a class of so-called isorefractive media, where the 
relation EeSReff = n? = EmRm holds. Thus, for effective 
permittivity we have 


e,^ = em[l-^y\ (4) 

Thus, we argue that for the correct description of the 
wire medium, Em and /im in Eqs. m should substituted 
with right hand-sides of Eq. (jT]) and Eq. (|3]), respectively. 
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We note, that there are a lot of works dedicated to the ar¬ 
tificial diamagnetism in metamaterialsi^^— . At the same 
time, surprisingly, diamagnetism in wire media metama¬ 
terials is conventionally abandoned. 

We also support the analytical expressions with the 
numerical modelling. For that we use the homogeniza¬ 
tion technique similar to the one described inii. We 
again assume the quasistatic approximation and consider 
the propagation of the TEM mode along the wires. In 
this case, we can set that electric field is aligned along 
y and magnetic field along x directions. We introduce 
static electric and magnetic potentials Ue, Um defined 
as = E(H). Then, two Laplace equations for 

Ue and Um can be solved. For the electric potential we 
use following boundary conditions: Ue\y=±d /2 = Ameo, 
due/dx\^^±d /2 = 0: and Ue\r=ro = 0 at the wire in¬ 
terface. For magnetic potential, boundary conditions 
are: Um\x=±d /2 = ±MmO, dumldm\y^±d /2 = 0 and 
dum/dr\r=ro = 0. We solve the both Laplace equa¬ 
tions numericallyi^. Knowing the magnetic potentials 
and therefore the electric and magnetic fieds, it is pos¬ 
sible to express the effective permittivity and magnetic 
permeability as 

I-d/2^">-^v\v=-d/'^^^ 
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The maps of the scalar electric and magnetic potentials 
as well as field lines are shown in Figl2Ia,b). We com¬ 
pare the analytical expressions Eqs. (I3]),(l4]) with the more 
rigorous approach given by expressions Eqs. (0,® in 
Figs.®c,d), respectively. We use the following structure 
parameters: £m = 16,/Tm = 1- It can be seen that for 
thin wires, the diamagnetism is negligible and both ap¬ 
proaches almost perfectly match. At the higher values 
of ro/d the two methods start to give slightly different 
results but still match quite well. The discrepancies at 
large ro/d appear, since the analytical approach consid¬ 
ers a single wire and does not account for the structure 
periodicity and interaction between wires in different unit 
cells, which grows with the growth of ro/d. 

Moreover, in order to check if the quasistatic ap¬ 
proximation is applicable in the case of realistic wire 
media structures, we have extracted the effective elec¬ 
tric and magnetic susceptibility using Nicholas-Ross- 
Weir Methodi^Ti^. We have calculated the reflection 
and transmission spectra of the wire medium slab with 
the same material parameters, thickness D = lOd in 
the full-wave electromagnetic simulation package^^.The 
wire radius rp was a parameter which was changed in 
the region [O.Old, 0.4d]. The conventional Nicolson-Ross- 
Weir (NRW) method essentially is an inversion of the 
Rayleigh expressions for the transmission and reflection 
from the uniform material slab with respect to effective 




FIG. 2. (a,b) Maps of the scalar electric (a) and magnetic (b) 
potentials and field lines for the electric (a) and magnetic (b) 
fields in the unit cell for the case of ro/d = 0.25, the color 
scale is the same for the both maps, (c) Effective dielectric 
permittivity vs ro/d obtained via the analytical formula in 
Eq. ([4|l (dashed black line), the semianalytical expression ([5]) 
(solid blue lines) and from extracting the effective parame¬ 
ters with NRW method (red circles), (d) Effective magnetic 
permeability vs ro/d obtained via the analytical formula Q 
(dashed black line), the semianalytical expression ® (solid 
blue lines) and from extracting the effective parameters with 
NRW method (red circles) 


material refractive index n = y/£eSldeS and impedance 
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The values of n, Z are then substituted to obtain 
the eeff,/reff- The results are shown with red dots in 
Figs. [2Kc,d). As shown in Figs. 0c, d) there is a good 
match between the full numerical simulation and semi¬ 
analytical approach. 

We stress that omitting the effective permeability of 
the structure leads to incorrect values of reflection and 
transmission coefficients if calculated with Rayleigh for¬ 
mulae, as shown in Fig. 0a), where the reflection and 
transmission spectra for the cases when the permeability 
is accounted for or omitted are compared to the numer¬ 
ical modelling. The value of ro/d is set to 0.25. We 
can see that omitting the magnetic response gives the 
incorrect results for the reflection and transmission spec¬ 
tra. Moreover, as can be seen in Fig. 0b), the extracted 
parameters correctly describe the electromagnetic prop¬ 
erties of the wire medium slab in the case of oblique in¬ 
cidence, when diamagnetism is taken into account. 
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FIG. 3. (a) reflection and transmission spectra for the case 
of normal incidence. The solid and dashed lines correspond 
to the reflection and transmission calculated with Eqs. 0,® 
with e = 16, p = 1 (dashed lines), and with e = 23.19,/r = 
0.69 (solid lines). The dots correspond to the CST modelling 
results, (b) The depedence of reflection and transmission on 
the incidence angle at the incident wave frequency 250 GHz; 
solid lines correspond to reflection and transmission calcu¬ 
lated with Eqs. 0,® with e = 23.19, = 0.69 and squares 

correspond to the CST modelling. Parameters of the wire 
medium structure are presented in the text. 


III. EXPERIMENTAL MEASUREMENTS 

In order to validate our results experimentally, we have 
constructed a wire metamaterial block with the follow¬ 
ing parameters: period of the structure d = 11 mm, wire 
radius rg = 5 mm, wire length D = 150 mm. The block 
consists of 20 x 20 wires with air dielectric matrix, i.e. 
Em = 1- The photo of the structure is shown in the 
inset of Figure S] The experimental measurements were 
conducted in a following manner: a rectangular horn an¬ 
tenna (TRIM 0.75 GHz to 18 GHz; DR) connected to 
a transmitting port of the vector network analyzer Agi¬ 
lent E8362C was used to approximate a plane-wave ex¬ 
citation. The metamaterial block was placed at the far- 
held region of the antenna and a similar horn antenna 
(TRIM 0.75 GHz to 18 GHz) was employed as a receiver. 
The effective scattering cross-section was obtained from 
the imaginary part of the forward scattering amplitude 
(due to the optical theorem). Due to the fact that the 
horn antenna modihes the wavefront substantially, we 
have normalized the scattering cross-section to unity. Af¬ 
ter completing the experimental measurements, we have 
performed the numerical simulation of the scattering at 
the metamaterial block of the same dimensions for three 
different models. In the flrst model, the metamaterial 
block was modelled as an array of perfectly conducting 
wires, reconstructing the geometry of the experimental 



FIG. 4. Scattering cross-section of the wire metamaterial 
block. The inset shows the photo of the structure. The ex¬ 
perimentally determined scattering cross-section is compared 
with numerical calculation performed in CST. The ripples in 
the experiment are explained in the main text. 


setup. In the second model, the block consisted of an 
uniform anistropic media with the following parameters: 
e_L = 5.762, £|| = -10000.0 -F l.Oi, = 0.1735, /ry = 1.0, 
where subscripts y and _l correspond to the directions 
parallel and perpendicular to the wire axis, respectively. 
The scattering cross-sections were normalized to unity. 
The above values for dielectric permittivities and mag¬ 
netic permeabilities have been derived with the expres¬ 
sions in Eqs. ©, ®. Finally in the third model, we have 
calculated the scattering on metamaterial block made 
from the anisotropic uniform media with the effective 
transverse dielectric and magnetic permeabilities equal to 
unity. Notably, it is quite obvious that in the third model, 
the scattering cross-section should be vanishingly small 
for the case of normal incidence, since both the effec¬ 
tive impedance and refractive index of this structure are 
equal to unity. The spectra of the forward scattering effi¬ 
ciency for the first two cases are shown in Fig. 2] The ob¬ 
served oscillations in the experimental scattering spectra 
are due to the artificial Fabri-Perot resonances between 
metamaterial block interface and source and receiving 
horn antennas. We observe a good correspondence be¬ 
tween the experimental cross-section spectra and numer¬ 
ical results, namely in the positions of the cross-section 
dips and in the cross-section modulation contrast. Even 
better correspondence is observed between the two nu¬ 
merical simulations. We note, that such correspondence 
could not be achieved in the approximation of effectively 
non-magnetic media, since as was discussed previously, 
such a media does not scatter the normally incident elec¬ 
tromagnetic radiation at all. The observed discrepancies 
between the experimental and numerical simulation re¬ 
sults are mainly due to the wavelength-dependent wave- 
front distortions of the incident wave, produced by the 
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horn antenna. 
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